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We introduce a new model of the discrete time quantum walk, the self-avoiding quantum walk, 
which is not allowed to step back onto a position which it has just occupied. This allows us to simulate 
a dimer and we achieve it by introducing a new type of coin operator. We describe its basic properties 
and provide numerical evidence that the standard deviation of the walker is constant regardless of 
the initial state. This contrasts strongly with previously studied coins such as the Grover operator, 
where the initial condition can be used to control the standard deviation of the walker. 



I. INTRODUCTION 

The discrete time quantum walk has been extensively 
studied since its introduction. Initially, various models 
were proposed HHS with a variety of applications in 
mind. In this paper we are concerned with the trans- 
port properties of new type of quantum walk. Quantum 
walks have been shown to have interesting transport 
properties in a variety of scenarios. On the line they 
achieve ballistic transport |4] and they were first shown 
to have an exponential speedup over the classical ran- 
dom walk on the hypercube by Kempe [5]. The faster 
hitting times have been used in algorithmic applications 

Ed. 

The quantum walks introduced thus far concern ide- 
alised walkers with no spatial extension. Whilst these 
have many uses in modelling physical and biological 
processes [8], in some cases we may want to consider 
walkers which do have spatial extension, and hence 
cannot move into space which they are already occu- 
pying. Self-avoiding random walks were developed to 
model precisely such processes, initially the folding of 
polymers. In some cases, such a restricted model may 
not be suitable. Continuing the example, polymers can 
overlap to form loops, but they cannot occupy volume 
that they already occupy. A type of intermediary be- 
tween the fully random walk and the self-avoiding walk 
is known as the non-reversal walk. A non-reversing ran- 
dom walker is allowed to occupy spaces on a lattice that 
it has previously occupied, but is not allowed to return 
to the point which it has just come from. In this paper we 
introduce a quantum version of such a walk, which ef- 
fectively simulates a dimer. The motivation for studying 
such a walk is much the same as that for studying the 
classical version: more realistic simulation of physical 
systems. 

In both the classical and quantum case, the self- 
avoiding or non-reversal walk on the line is trivial. This 
is because there are only two degrees of freedom in the 
movement, so if one of those is prohibited by the model, 
then unidirectional ballistic transport is obtained. The 
walks studied in this paper take place over a square lat- 
tice, in which case the dynamics are highly non-trivial. 



The paper proceeds as follows: In Section I A classical 
self-avoiding walks are introduced in more detail. Then 
for the sake of comparison the properties of the quan 
turn walk on the square lattice are discussed. The self- 
avoiding quantum walk is then defined in Section II A 



and its properties are described in Section II B 



A. Classical self-avoiding random walks 

The classical self-avoiding walk has proven difficult to 
treat analytically, hence the results concerning it have all 
so far been numerical [9 J and there remain many open 
questions about it. Even enumerating the number of self- 
avoiding walks has proven very difficult, despite them 
being so rare that coming upon one by mistake when 
examining a random walk is highly improbable. If c n is 
the number of walks of precisely n steps, then the total 
number of self-avoiding walks up to length n is Yun>2 c n- 
Some facts are clear, for example that c n+m < c m + c n . 
The set of walks of length n concatenated with the set 
of walks of length m contains not only the self -avoiding 
walks of length n + m but some which overlap, hence the 
inequality. However, determining the precise number of 
walks is difficult, though bounds have been established. 
The number of self-avoiding walks must be less than 
the number of non-reversal walks, as these include the 
self-avoiding walk as a subset. Additionally it is pos- 
sible to construct subsets of self-avoiding walks which 
grow as 2 n , so we know that there are between 2 n and 
3 n self -a voiding random walks. The best evidence so 
far suggests that the number of self-avoiding walks of 
length n is proportional to 2.638 w , and this is provided 
as a non-rigorous estimate in [10]. The best evidence for 
this value was obtained by enumerating each such walk 
of length up to 51 and required a 1024 processor super- 
computer lUTl . Without new algorithms it is unlikely 
that we will be able to enumerate much further than 
this. Guttmann obtained 14,059,415,980,606,050,644,844 
walks using this method, however these walks make up 
less than 1 in 240 million of the possible random walks 
of this length. 

As even counting the walks has proved difficult, it is 
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FIG. 1. (Color online) Probability distributions arising from a) 
the DFT coin and b) the Grover coin for the symmetric initial 
state shown in equation|S] 



unsurprising that little is known regarding other prop- 
erties. The property we are most interested in when 
comparing walks is the standard deviation. The mean 
squared displacement is conjectured to be n 3/2 though so 
far, even a proof that the exponent must be between 1 
and 2 is elusive l9l . 

Another known fact about self-avoiding walks 
demonstrates a key difference between the self-avoiding 
walk and its standard and non-reversal counterparts. 
This is that the self -avoiding walk does not necessarily 
continue indefinitely. This is because it is possible to 
reach a lattice site whose only adjacent lattice sites have 
previously been visited, hence the walker cannot con- 
tinue. On any given step there is a small probability that 
the walk will end. 

The non-reversal walk is in some ways more tractable. 
For example, it is clear that on the square lattice there are 
3 n such walks, where n is the number of steps traversed 
by the walker. Its mean squared displacement is 2n, so it 
spreads twice as fast as the standard random walk, but 
slower than the completely self-avoiding walk. There 
is very little literature on the non-reversal walk, and 
this tends to examine specific characteristics of the walk 
relevant to the study of polymer chains [12J, rather than 
its general features. 



B. Quantum walks on the square lattice 

The properties of the discrete time quantum walk on 
the square lattice were extensively explored in [ 13] fol- 
lowing initial investigations in tl4l . In particular they 
examined the mean position at time t: 

(r) t = Y J ylx 2 + y 2 p(x,y,t) (i) 

x,y 

so r is the radial distance from the origin, and y, t) is 
the probability of finding the walker at position (x, y) at 
time t. They also characterise the walks in terms of the 
standard deviation: 

a = V<r 2 ) - (r) 2 (2) 

The authors use three choices of coin operator. Their first 
choice behaves like the Hadamard operator for both the 
Teft/right' (|/)/|r)) coin states and the 'up/down' (\u)/\d)) 
coin states. This operator does not mix the two dimen- 
sions, so a two dimensional version of the distribution 
of the walk on the line is obtained. More interestingly, 
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FIG. 2. (Color online) Probability distributions arising after 100 
steps of a typical self-avoiding quantum walk shown a) over 
the entire lattice and b) in profilep] 
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they consider the DFT and Grover operators for a num- 
ber of initial conditions. These coins are both unbiased, 
in that they distribute amplitude equally between each 
coin state. In contrast to the walk on the line, they find 
that the dynamics of the walk depend strongly on the 
coin used. Additionally, the dynamics for a specific coin 
depend strongly on the choice of initial state. The lowest 
and highest standard deviations obtained for the posi- 
tion of the walker were found using the Grover operator. 
It was observed that the reason for this is that regardless 
of the initial state, the distribution forms a central spike, 
with a ring around it which propagates outwards. The 
choice of initial condition controls how much amplitude 
is situated in the central spike, and how much amplitude 
is situated in the ring. The initial state at the origin for 
both the distributions plotted in Figure [T] is: 
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(3) 



The authors of 1 13] studied all unbiased four dimen- 
sional unitary operators with entries equal to either +1/2 
or ±//2, which when the leading diagonal is selected to 
be 1/2 gives 640 unitary operators. These operators were 
found to fall into ten types, with the DFT, Hadamard and 
Grover all being of different type. Of all operators tested, 
the Grover operator was found to give rise to both the 
largest and smallest standard deviation, depending on 
the initial state. 



II. SELF- AVOIDING QUANTUM WALKS 

A. Definition 

Our definition of a self-avoiding quantum walk is 
more similar to the classical non-reversal walk than the 
fully self-avoiding walk. We require that the diagonals 
of the coin operator are all equal to zero, so that ampli- 
tude which has come from a given node cannot return 
to that node during the next step of the walk. A general 
coin fulfilling this criteria is specified by: 
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where all of the variables are real and < y 2 + A 2 < 1, 
and /(A, y) = ^1 - (A 2 + y 2 ). A simple example of such a 
coin, used to produce the probability distribution shown 
in Figure pi takes 6 = <p = n and A = y = f(A,y) - 4=, 

with all other phase variables being zero. This leads to 
the following coin: 
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so the initial amplitude in each coin state is redistributed 
between the other coin states, and none of it is transmit- 
ted to the node it came from. 

Whilst the walks examined below use real parameters, 
it is clearly possible to have coins with complex numbers 
in them, such as 



We can see that this coin does what it is required to by 
examining its operation on an arbitrary state 
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It is clear that 



but the properties discussed below have been verified 
numerically for coins with real entries only. 

A self-avoiding walker of length three, which models 
a trimer, requires a more complex definition. This is be- 
cause the trimer on a square lattice, presuming that the 
head and tail are distinguishable, can be in one of three 
configurations, as shown in Figure [3] These configura- 
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FIG. 3. 
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The three possible configurations of a self avoiding 



tions can then have four possible orientations, facing in 
the up, down, left and right directions. This means that 
in order to represent the trimer, a qutrit tensored with a 
ququat is required, giving rise to a twelve dimensional 
coin. The head of the trimer can make the same moves 
as that of the dimer, but now as well as updating the 
position of the walker, the configuration, which will in 
general be in a superposition, is also updated. 

Whilst it may be possible to define a fully general self- 
avoiding quantum walk, this will in practise be difficult. 
Even deciding when a site has been visited has some 
ambiguity, for instance, can one part of the superposition 
access sites which another part of the superposition has 
previously accessed? However, the self-avoiding coin 
developed here has sufficiently interesting properties for 
it to be worth exploring regardless, and we now turn to 
these. 



B. Properties 



An example of a typical probability distribution aris- 
ing from a self-avoiding quantum walk is shown in Fig- 
ure [2] The distribution is shown in 3d and in profile 
in order to show the similarity of the profile to that of 
the walk on the line. Regardless of initial condition, the 
dynamics are similar, tracing out a square with peaks at 
each corner. The initial condition determines the height 
and number of distinctive peaks at the corners of the 
square. The initial condition used to obtain the distri- 
bution shown in Figure [2] is the same as that used in 
Equation [3] 

These walks were investigated by varying the param- 
eters in the coin as well as the initial condition. It was 
found numerically that for all choices of initial condition 
the mean position (Equation [l| and standard deviation 
(Equation|2l of the walker is constant at a given time t, as 
can be seen for the standard deviation in Figure [4] Fur- 
ther investigations suggest that all the moments of the 
distribution are independent of the initial condition. The 
parameters which determine the standard deviation are 
A and y, so these values determine some real constant 
k for each self-avoiding coin so that the standard devi- 
ation at time t is equal to kt. This contrasts strongly 
with standardly used coins, as mentioned in Section I B 
above, where the initial condition can be used to control 
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FIG. 4. Comparison of standard deviation for an arbitrary 
initial state for (a) the Grover coin and (b) the self-avoiding 
coin. 



the standard deviation of the walk using the Grover coin 
operator. The mean and standard deviations as a func- 
tion of time are shown for a variety of coins in Figure [5] 
Whilst the mean position increases faster than all other 
coins apart from the Grover, the standard deviation in- 
creases more slowly, so in cases where we have freedom 
to choose our initial condition, the Grover operator is still 
the best to use if we want to achieve the fastest possible 
spreading of the amplitude across the lattice. 



III. CONCLUSION 

We have introduced a new type of coin for the dis- 
crete time quantum walk and shown numerically that 
it has some notable properties, namely that the mean 
and standard deviation of the position of the walker is 
independent of the choice of initial condition. Ideally, 
an analytic proof that these properties do not depend 
on the initial condition is required, as well as an ana- 
lytic expression for the standard deviation in terms of 
the relevant coin parameters, but these have eluded us 
thus far. The standard deviation grows linearly with t 
for much the same reason as it does the walk on the line, 
as the coin operator always ensures that some amplitude 
moves away from the origin with each step. However, 
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FIG. 5. (Color online) Comparison of coin operators in terms 
of a) mean and b) standard deviation. The initial condition 
leading to the steepest gradient for each coin was used. 



this does not explain why the standard deviation is de- 
termined by the choice of A and y alone. 



A fully self-avoiding quantum walk is yet to be de- 
veloped and explored, and thus far no continuous time 
analogues have been proposed. In order to fully capture 
the folding of polymer chains, the quantum nature of 
their constituents may need to be accounted for, which 
would provide a potential application for such a model. 

In addition to further investigating the properties of 
the self-avoiding quantum walk and developing the con- 
cept, better comparisons with its classical counterpart 
may be useful. For example, only walks on square lat- 
tices have been considered here, but the self-avoiding 
random walk has been shown to have macroscopic prop- 
erties which are independent of the choice of lattice. In 
order to see if this property carries over into the quantum 
case, coins of varying dimension are required. 

In short, there are many potential avenues through 
which to explore the self-avoiding quantum walk, and 
the choices of which direction to take will depend on 
the desired applications of the walk. The purpose of the 
work here was simply to find out whether using such a 
self-avoiding coin operator gives rise to any properties 
which are not observed using other operators, hopefully 
provoking further investigation into the uses of such a 
coin. 

ACKNOWLEDGMENTS 

KB is funded by the UK Engineering and Physical Sci- 
ences Research Council. VK is funded by a Royal Society 
Research Fellowship. TP, BH, and AB received support 
from Royal Society summer bursaries and EPSRC.VP 
was supported by IASTE and SM by JFT Grant ID 15619 
and ANR-10-JCJC-0208 CausaQ grant. . 



[1] S. Gudder, Quantum Probability. Academic Press Inc., CA, 
USA, 1988. 

[2] Aharonov, Y„ Davidovich, L., and Zagury, N., "Quantum 
random walks," Phys. Rev. A, vol. 2, p. 16871690, 1992. 

[3] G. Grossing and A. Zeilinger, "Quantum cellular au- 
tomata," Complex Systems, vol. 2, p. 197 208, 1988. 

[4] A. Ambainis, E. Bach, A. Nayak, A. Vishwanath, and J. 
Vishwanath, "One-dimensional quantum walks," Proc. 
33rd Annual ACM STOC, pp. 60-69, 2001. 

[5] J. Kempe, "Quantum random walks hit exponentially 
faster," Probability Th. and Related Fields, vol. 133 (2), 
p. 215235., 2005. 

[6] A. Childs, R. Cleve, E. Deotto, E. Farhi, S. Gutmann, and D. 
Spielman, "Exponential algorithmic speedup by quantum 
walk," Proc. 35th annual ACM STOC, pp. 59-68, 2003. 

[7] A. Childs, E. Farhi, and S. Gutmann, "An example of the 
difference between quantum and classical random walks," 
Quantum Inf. Process., vol. 35, no. 1, 2002. 

[8] M. Mohseni, P. Rebentrost, S. Lloyd, and A. Aspuru- 
Guzik, "Environment-assisted quantum walks in pho- 



tosynthetic energy transfer," /. Chem. Phys., no. 129, 
p. 174106, 2008. 
[9] B. Hayes, "How to avoid yourself," AmSci, vol. 86, p. 314, 
1998. 

[10] I. Jensen, "A parallel algorithm for the enumeration of 
self-avoiding polygons on the square lattice," /. Phys. A: 
Math. Gen, vol. 36, pp. 5731-5745, 2003. 

[11] A. J. Guttmann and A. R. Conway, "Square lattice 
self-avoiding walks and polygons," Ann. Comb., vol. 5, 
pp. 319-345, 2001. 

[12] A. Skrilos and S. Chirikjian, "Position and orientation 
distributions for non-reversal random walks using space- 
group fourier transforms," / Algebr Stat., vol. 1, pp. 27-46, 
2010. 

[13] B. Tregenna, W. Flanagan, R. Maile, and V. Kendon, "Con- 
trolling discrete quantum walks: coins and initial states," 
New J. Phys., vol. 5, no. 83, 2003. 

[14] T. D. Mackay, S. D. Bartlett, L. T. Stephenson, and B. C. 
Sanders, "Quantum walks in higher dimensions," /. Phys. 
A: Math. Gen., vol. 35, 2002. 



